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EFFECTS OF PRECESSION AND NUTATION 


by 

James P. Murphy 


SUMMARY 


Transformations to account for the precession and nutation in the equator 
and equinox of the earth are given in units of degrees with time measured in 
Julian days and initial epoch of 2436099.5 J.D., the Julian Date for Space. 
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EFFECTS OF PRECESSION AND NUTATION 


by 

James P. Murphy 


INTRODUCTION 

The true equator and equinox of the earth are in continuous motion in space. 
The motion consists of two parts, general precession and nutation. General 
precession carries the mean pole of the equator around the pole of the ecliptic 
in about 26,000 years and nutation carries the true pole around the mean pole in 
about 18.6 years. Thus the effect of precession is secular and the effect of nuta- 
tion is periodic. 

Three separate coordinate systems are to be considered. The Z -system is 

r \j 

moving under the influence of precession and nutation; the Z -system is moving 
under the influence of nutation; and the W-system is fixed in space. These co- 
ordinate systems are now defined. 

hi the Z -system the Z 3 axis points toward the instantaneous direction with 
respect to the celestial sphere of the mean polar axis, and Z points toward the 
instantaneous vernal equinox. At the same time in the Z -system, the Z 3 axis 
points toward the mean direction of the mean polar axis, and Z points toward 
the mean vernal equinox. At some specified epoch in the W-system, the W 3 axis 
points toward the mean direction of the mean polar axis and W x points toward 
the mean vernal equinox. 

PRECESSION 

Let k and u> be the angle between the line of intersection of the mean 
equator at T c and the mean equator at T with the W 2 and lt 2 axes, respectively 
(See Figure 1). Let v be the angle between the planes of the two mean equators. 

The relationship between the Z and W coordinate systems then is 

Z = R 3 (- to) R 2 (v) R 3 (-k) W 

and (1) 

W = R 3 (k) R 2 (- v ) R 3 0) Z 
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Figure 1 


where R (a) indicates a rotation matrix with the rotation about the n th axis by 
an angle a. 

The explicit relationships indicated by Equations (1) are 
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The values of k, co, and v are 

k = [23042 "53 + 139'.'73 (t o - 1.900) + 01*06 (T o - 1.900) 2 ] (T - Tj 

+ [301*23 - 01*27 (T o - 1.900)] (T - T c ) 2 + 181'00 (T - T o ) 3 

co = [230421*53 + 1391*73 (T o - 1.900) + 01*06 (T o - 1.900) 2 ] (T - T o ) 

+ [l091'50 + 01*39 (T o - 1.900)] (T - T o ) 2 + 181*32 (t - T o ) 3 

v = [200461*85 - 851*33 (T o - 1.900) - 0*.'37 (t - 1.900) 2 ] (T - T o ) 

+ [- 421'67 - 01'37 (T o - 1.900)] (t - T o ) 2 - 41*180 (T - T o ) 3 

where T 0 and T are measured in thousand tropical years and where T 0 refers 
to the epoch time of the fixed mean equator and T refers to the time of the 
moving equator. 

NUTATION 

Let A/x, t\v , and Ae be the nutation in right ascension, declination, and 
obliquity, respectively. The relationship between the Z and Z systems is then 

Z = R 3 (- A/x) R 2 (Av) Rj(- Ae) Z 

and (2) 

Z = R x (Ae) R 2 (- Av) R 3 (A/x) Z 

The explicit relationships indicated by Equations (2) are then 


z, 

z 2 


cos A V cos A y - cos A (i sin A v sin A e - cos A e sin A 

cos A V sin A y. - sin A y. sin A v sin A e + cos A e cos A y 

sin A v sin A e cos A v 


- cos A y sin A v cos A e + sin A e sin 

- sin A y sin A v cos Ae - sin A e cos 
cos A c cos A v 



z i 




z 2 


'V 


_ Z J 
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and 


m 

z = 

L z iJ 

The expressions for Aft, Av, and A e may be obtained from expressions for 
the nutation in longitude A \p and in obliquity A e appearing on pages 44 and 45 of 
Reference 3. 

This relationship is illustrated in Figure 2. 


cos A v cos A (j. cos A p sin A /j, sin A v 

' cos A y sin A v sin A e - cos A e sin A y - sin A y sin A v sin A e + cos A e cos A y sin A e cos A v 

• cos A y sin A v cos A e + sin A e sin A y - sin A /x sin A v cos A e - sin A e cos A y cos A € cos A v 


jlAj 



Figure 2 


Thus, 


A/i. - Ai/> cos e 
Av = A</< sin e 


(3) 


The "equation of the equinoxes" or nutation in right ascension is the right 
ascension of the mean equinox referred to the true equator and equinox and 
represents the difference between the mean and true right ascensions for a body 
on the equator. It is thus the difference between mean and apparent sidereal 
time. 


4 


FORMULAS FOR EPOCH OF THE JULIAN DATE FOR SPACE 


t 


If we wish to use the Julian date for space as the initial epoch, the day as the 
time unit and the degree as the unit of angular measure, then some formulas of 
the previous sections undergo changes in the coefficients. That is, the preces- 
sion formulas for k, co, and v are adjusted for the new units. An adjustment 
must be made in the coefficients in the series for Ai p and Ae as well as the 
motion of the arguments of the various terms. 

The formulas for precession are then given by 

k = [l?753067 x IO' 5 + 2?9097 x IO -13 (t Q - tj D s ) 

+ 3?4 x IO' 22 (t o - t J-D . s .) 2 ] (t - tj 
+ [fi?291 x ur 14 - 1?5 x IO' 21 (t o - tj. D>s> )] (t - tj 2 

+ 1?026 x 10" 19 (t - t o ) 3 

a = [l?753067 x 10" 5 + 2?9097 x 10' 13 (t o - t ; D s ) 

+ 3.4 x 10 22 (t Q - tj D s ) 2 J [t ~ t 0 ) 

+ [2?2805 x 10" 13 + 2?3 x 10~ 21 (t Q - tj D s )] (t - tj 2 

+ 1?044 x 10‘ 19 (t - t o ) 3 

v = [ 1.524289 x 10' 5 - 1?7777 x 10 -13 (t o - t J-D#s J 

- 2?1 x 10- J1 (t„ - t J D S .) 2 ] (t-tj 

+ [- 8?889 x 10‘ 14 - 2fl x IO’ 21 (t 0 - t J D S .)] (t - t a 2 ) 

- 2? 383 x IO" 19 (t - tj 3 

where tj. D .s. = 2436099.5 J.D. and where (t D - tj. D .s.) and (t - t D ) are 
measured in Julian days. 

In the series for nutation in longitude and obliquity terms of the order of a 
hundredth of a second of arc or larger are kept. This amounts to about one 
quarter of the terms listed in Reference 3 for these series. Thus 
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A«/> = 


(- 4.78965 x 10“ 3 - 1.323 x 10“ 10 t) sin Q 
+ (5.800 + 10~ 5 + 1. x 10' 13 t) sin 2Q 

+ (- 3.5361 x 10“ 4 - 1. x 10" 12 t) sin 2 (F - D + 0) 

+ (3.497 x 10" s - 2.4 x 10~ 12 t) sin £’ 

+ (- 1.38 x 10~ 5 + 9. x 10“ 13 t) sin ( £' + 2F - 2D + 20 . ) 

+ (- 5.94 x 10“ 6 + 4. x 10‘ 13 t) sin (£* - 2F + 2D - 2D) 
+ (3.44 x 10" 6 + 8. x 10“ 14 t) sin (2F - 2D + Q) 

+ (- 5.659 x 10“ 5 - 1. x 10“ 13 t) sin 2 (F + 0) 

+ (1.88 x 10" s + 8. x 10“ 14 t) sin £ 

+ (- 9.50 x 10" 6 - 3. x 10“ 13 t) sin (2F + ft) 

- 7.25 x 10“ 6 sin (£' + 2F + 2 Cl) 

- 4.14 x 10“ 6 sin (£'- 2D) 

- 3.17 x 10“ 6 sin (£'- 2F - 2fi) 

+ 1.7 x 10 -6 sin 2D 

+ 1.6 x 10” 6 sin (£' + fl) 

+ 1.6 x 10“ 6 sin (£’ - f2) 

+ 1.4 x 10 -6 sin (£' - 2F - 2D - 2fi) 
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and 

Ae 


(2.5585 x 1C)' 3 + 6.9 x 10" 12 t) cos 0 

+ (- 2.51 x 10" 5 + 3. x 10“ 13 t) cos 2 0 

+ (1.533 x 10' 4 - 2.2 x 10” 12 t) cos 2 (F - D + 0 ) 

+ (5.99 x 10 -6 - 4.4 x 10" 13 t) cos (£' + 2F - 2D + 2D) 

+ (- 2.58 x 10~ 6 + 2. x 10” 13 t) cos (£ - 2F + 2D - 2 0) 

- 1.8 x 10" 6 cos (2F - 2D + 0) 

+ (2.46 x 10” 5 - 3.9 x 10" 13 t) cos 2 (F + 0) 

+ 5.08 x 10” 6 cos (2F + Q) 


+ (3. 14 x 10” 6 - 8. x 10” 14 t) cos (£ + 2F + 20) 

- 1.4 x 10” 6 cos ( l - 2F - 20) 

where 

l - 299?615949 + 13?0649927374d + 6?909 x 10 _12 d 2 + 2?95 x 10“ 19 d 3 

V = 254?436609 + 0?9856002621d - 0.116 x 10“ 12 d 2 - 0f68 x 10” 19 d 3 
F = 239?342609 + 13?2293503475d - 2?407 x 10“ 12 d 2 - 0?07 x 10” 19 d 3 
D = 285?634588 + 12? 1907491461 d - 1?074 x 10" 12 d 2 + 0?39 x 10" 19 d 3 
0 = 222.941764 - 0?0529538565 d + 1?560 x 10' 12 d 2 + 0?46 x 10' 19 d 3 


( 5 ) 


and where d is measured in Julian Ephemeris Days since the Epoch Date of the 
Julian Day for Space, 2436099.5 J.D. 
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APPENDIX A 


LIST OF SYMBOLS 


Z- system 


Z-system 

W-system 

K, 00, V 

A/x 

Av 

Ae 


A 0 

T 

T' 


T 0 

T 


t 

^ J . D . S . 

R„(“) 

i, F, D, n 


coordinate system moving under the influence of precession 
and nutation 

coordinate system moving under the influence of nutation 

coordinate system fixed in space 

precession angles given by Equations (4) 

nutation in right ascension 

nutation in declination 

nutation in obliquity 

nutation in longitude 

true equinox 

mean equinox 

epoch time of fixed mean equator in tropical centuries 

epoch time of moving equator in tropical centuries 

T o expressed in Julian days 

T expressed in Julian days 

Julian Date for Space = 2436099.5 J.D. 

rotation about the i * h axis by an angle a 

angular variables defined in Equations (5) 


A-l 



